1 


Comment on "Multiferroicity Induced by Dislo- 
cated Spin-Density Waves" 

A recent Letter by Betouras et al. jl| (BGB) proposed 
a phenomenological magnetoelectric (ME) interaction to 
explain the spontaneous polarization P in a class of ma- 
terials of which YMn2 05 (YMO) was taken to be the best 
exemplar. They claimed 1) its behavior agreed perfectly 
with their theory and 2) previous theories failed in this 
regard. Here I show that these claims are incorrect and 
that their interaction is unimportant in most systems. 

Figure 1 shows P measured by Kagomiya et al. J3| 
compared to the microscopic theory of Chapon et al. [3] ■ 
The ME interaction proposed by BGB gives P=0 in the 
low-temperature (T < 20K) incommensurate (IC) phase 
and, as I will discuss, nearly zero in the high-temperature 
(T > 20K) commensurate (CM) phase. Below I show 
that the Landau approach 5] fits the data qualitatively. 


Fig. 1. (Based on Ref. 
i.) P vs T for YMO. 
Filled circles: theory of 
Ref. i. Thin line: data 
of Kagomiya et al. The 
dashed line indicates the 
IC-CM phase boundary. 
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The results of Ref. 1 follow from a gradient expan- 
sion of the free energy to mimic the variation of mag- 
netization on an atomic length scale. However, the most 
successful long- wavelength treatment of antiferromagnets 
[3| utilizes the symmetry of both the staggered and direct 
magnetization. The generalization to the multisublattice 
case of YMO is to consider all the symmetry elements of 
the unit cell as done in Refs. Si. There one introduced 
order parameters (t„ of known symmetry to describe the 
magnetism of the unit cell and for YMO the trilinear 
ME coupling is 3 F = aP[|crip - |cr2p], which gives 
P — const X [|(T2p - IcTiP]- I will caU this the "standard" 
model and it was used Q to successfully interpret the 
polarization of the CM phase of YMO. I show here that 
this formulation can also explain the smaller nonzero po- 
larization in the IC phase which contradicts Ref. [H. For 
this purpose we need to extract the orderparameters for 
YMO since they were not given in Ref. |3- The symme- 
try of the spin pattern shown in their Fig. 2a [iQ] for 
the CM phase indicates that the order parameter called 
(Ti m Ref. i is zero. [9] To normalize the results, I de- 
fine fT2 = 1. Although the size of the moments does not 
change much at the CM-IC transition, the orientations 
of the moments do change and this change induces a sig- 


nificant redistribution of the order parameters. To see 
this consider the relation in their Fig. 2a between the 
Mn3+ spins at (0.09,-0.15,0.5) and at (0.41,0.35,0.5). In 
the CM phase their si-components are the same but their 
y-components change sign. In the IC phase both x and 
y- components change sign. This change in spin orienta- 
tion induces significant changes in the order parameters. 
I found the order parameters of the T — 1.9K structure 
to be such that |(T2p — Wi]"^ lies between -0.25 and -1-0.25. 
Even with this uncertainty it is clear that the change in 
the magnetic structure can explain the precipitate de- 
crease in P at the CM to IC phase boundary. 

Since the microscopic model of Ref. i is so satisfactory, 
one can ask whether it is consistent with either or both 
of the competing phenomenological theories. To answer 
this question I write the ME interaction of BGB in the 
language of order parameters effectively as 


V = const X PaiqY'Sinq, G) 


(1) 


where the S function conserves momentum to within a 
nonzero reciprocal lattice vector G (so that nq = G). 
For YMO, ri = 4 in Eq. (1), since one of the magnetic 
unit cell axes is four times that of the paramagnetic unit 
cell. The theory of Chapon et al. is quadratic in the 
spin operators. It can not give rise to 4q — G, because 
this would require either a four-spin ME interaction or 
treating the standard two-spin ME interaction in second- 
order perturbation theory. However, this latter possibil- 
ity is quite remote because the ME interaction is already 
quite weak in systems where the ME interaction acts in 
first order perturbation theory. So the Chapon et al. 
model, which fits the data nicely, is not a microscopic 
representation of the ME interaction of BGB. But since 
the model of Ref. 2 gives the polarization as a quadratic 
form in the spin operators, it is the microscopic analog of 
the "standard" model, which I therefore claim does rep- 
resent the data well. Moreover, rather than happening in 
a wide class of materials, the BGB interaction might be 
important, but only when the magnetic unit cell involves 
doubling of the paramagnetic unit cell. 

In summary: the existing "standard" model explains 
the data of Fig. 1 qualitatively, whereas that of BGB 
does not and is unlikely to be important for most systems. 

I thank A. P. Ramirez and T. Yildirim for suggestions. 
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° The top (bottom) panel should be labeled 24.7K (1.9K). 


